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ABSTRACT 

The bright Seyfert 1 galaxy MCG-6-30-15 shows large variability on a variety of time scales. We study the 
^ 3 day time scale variability using a set of simultaneous archival observations that were obtained from RXTE and 
the Advanced Satellite for Cosmology and Astrophysics (ASCA). The RXTE observations span nearly 10 6 sec and 
indicate that the X-ray Fourier Power Spectral Density has an rms variability of 16%, is flat from approximately 
10~ 6 -10~ 5 Hz, and then steepens into a power law oc f~ a with a <L 1 . A further steepening to a sa 2 occurs between 
10 -4 — 10~ 3 Hz. The shape and rms amplitude are comparable to what has been observed in NGC 5548 and Cyg X- 
1, albeit with break frequencies that differ by a factor of 10~ 2 and 10 4 , respectively. If the break frequencies are 
indicative of the central black hole mass, then this mass may be as low as 10 6 M . An upper limit of ~ 2 ks for 
the relative lag between the 0.5-2 keV ASCA band compared to the 8-15 keV RXTE band was also found. Again 
by analogy with NGC 5548 and Cyg X-l, this limit is consistent with a relatively low central black hole mass. 

Subject headings: galaxies: individual (MCG-6-30-15) — galaxies: Seyfert — X-rays: galaxies 
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1. INTRODUCTION 

The type 1 Seyfert galaxy MCG-6-30-15 has in recent years 
been the subject of intense study owing to the discovery by the 
Advanced Satellite for Cosmology and Astrophysics (ASCA) of 
a resolved , broad iron Ka flu orescent line in its hard X-ray 
spectrum (Ianaka et al. 1995). The shape of the line is con- 



mains constant. Iwasawa et al. (1999) suggest that multiple 



sistent with a gravitationally and Doppler shifted emission line 
which originates from near the inner edge of an accretion disk 
around a black hole. In the X-rays, MCG-6-30-15 is also one 
of the brighter and more variable type 1 Seyferts. It is therefore 
hoped that by examining the correlated variability between the 
ionizing X-ray continuum and various components of the Fe 
Ka line, one can obtain a size scale for the sys tem and thence 
a mass for the black hole (Reynolds et al. 1999). This requires, 
however, that the various components of the line be spectrally 
resolved on time scales shorter than the intrinsic response time 
scale of the line-emitting material. 

At present, ASCA, BeppoSAX, and Chandra are the only X- 
ray telescopes with the spectral resolution to measure fluxes in 
different line components separately. Unfortunately, the inte- 
gration times (;> lOks) required to obtain this resolution are 
longer than the light travel time across the inner edge of an 
accretion disk around a 1O 8 M Schwarzschild black hole. Fur- 
thermore, the fits to t he iron line appear to require a Kerr ge- 
ometry (although see Reynolds & Begelman 1997) in order to 
explain the broad red wing of the line and the lack of signif- 
icant emission blueward of 6.4 keV, thus making the relevant 
time scales even shorter. Nonetheless, over longer time scales, 
ASCA has measured signifi cant variability in th e shape of the 
MCG-6-30-15 Fe Ka line ( [Iwasawa et al 1996|). 

Time-resolved spe ctral investigations ( [iwasawa et al. 1996 ; 
Iwasawa et al. 1999|) suggest that the broad and narrow com 



ponents of the iron line are correlated differently with the flux 
state depending on the time scale investigated. For integrations 
> 10 4 s, the narrow component varies with the continuum flux 
whilst the broad component appears to be anti-correlated. In 
contrast, on shorter time scales the broad component responds 
immediately to flux changes whilst the narrow component re- 



localized X-ray flares occur on the disk surface near the inner 
edge and illuminate only a relatively small region of the disk. 
These small regions make contributions to narrow ranges in line 
redshift and thus produce very com plex temporal behavior . 



In support of this interpretation, Iwasawa et al. (1999) 



con- 
sider the iron line associated with a very short, bright flare 
from MCG-6-30-15 observed during 1997 August by ASCA 
(see Fig. U). The line was very redshifted and had little or no 
emission blueward of 6 keV. Its shape was consistent with hav- 
ing originated entirely from a small region at r ~ 5 GM/c 2 on 
the approaching side of the disk. The flare lasted about 4 ks, 
which is also approximately the orbital period at this radius for 
a 10 7 M Q black hole. In order for the line not to be significan tly 
more smeared by the orbital motion, Iwasawa et al. (1999) es- 
timate a black hole mass of ~ 2 x 1O S M . Alternatively, they 
suggest that a much lower mass is possible if all the line emis- 
sion arises from r < 5 GM/c 2 . 

A major complication that a high-mass model faces is the 
large amplitude and rapid X-ray variability of MCG-6-30-15 
(Fig. nh. In par ticular, during the perf ormance verification 
phase of ASCA, [Reynolds et al. ( 1 995)] noted that the 0.5- 
lOkeV flux increased by a factor 1.5 over 100s, i.e., the dy- 
namical time scale at the inner disk edge surrounding a max- 
imal Kerr, 2 x 10 6 M Q black hole. Typically one expects the 
bulk of the X-ray variability to occur on dynamical time scales 
or longer. It is clearly important to determine whether this vari- 
ability represented a rare, rapid event or if such time scales are 
truly characteristic of the behavior of MCG-6-30-15. 

One might expect that characteristic time scales should scale 
wi th the mass of t he c entral object (see, e.g., t he discussions 
of |McHardy 1988| and |Edelson & Nandra 1999| ); therefore, in 
this work we try to gauge the size of the system and the mass 
of the black hole in MCG-6-30-15 by studying its character- 
istic X-ray variability properties in comparison to other black 
hole systems such as Cyg X-l and NGC 5548. For the high- 
frequency variability analysis, we use the 1997 August simulta- 
neous ASCAI Rossi X-ray Timing Ex plorer (R XTE) observatio n 
discussed by Iwasawa et al. (1999) (see also Lee et al. 1998). 
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Fig. 1. — The RXTE 8-15 keV band (top) and the ASCA 0.5-2 keV band (bottom) lightcurves, in 4096 s timebins. These data have been obtained 
from the archival 1997 August simultaneous ASCAIRXTE observation (Lee et al. 1998). 



In our analysis we screen the ASCA data as outlined by Brandt 
et al. (1996), except that we use the more stringent criteria of 
7 GeV/c for the rigidity and an elevation angle of 10°. Data 
from both SIS detectors are combined into a single lightcurve. 
For the RXTE data, we use screening criteria and analysis tech- 
niques ap propriate for faint sou rces, as we have previously dis- 



cussed in phiang et al. (2000J , Specifically, we only analyze 



top layer data from proportional counter units 0, 1, and 2. 

2. POWER SPECTRA 

MCG-6-30-15 is detected by the All Sky Monitor (ASM; 



Levine et al. 1996) on RXTE with a mean countrate of 0.44 cps 
in the 1.3-12.2keV band. Currently, there is a 0.1 cps uncer- 
tainty in the zero level offset, as well as comparable magni- 
tude systematic time-dependent variations due to, for exam- 
ple, the solar angular position relative to the source (Remil- 
lard 1999, priv. comm.). Assessing the very low frequency 
(/ < 10" 6 Hz) X-ray variability of MCG-6-30-15 is therefore 
problematic. Efforts are currently underway, however, to re- 
vise the ASM data reduction process in order to minimize such 
systematic effects (Remillard 1999, priv. comm.); therefore, an 
ultra-low frequency variability study of MCG-6-30-15 in prin- 
ciple will be feasible in the near future. 

We are able, however, to investigate the high-frequency 
power spectral density (PSD) of MCG-6-30-15 by using the 
RXTE 8-15 keV and the ASCA 0.5-2 keV lightcurves binned 
on 4096 sec time scales. Both lightcurves contain data gaps, 
especially on the orbital time scale of ss 5 ks due to blockage 
by the Earth and passage through th e South Atlant ic Anomaly. 
We th erefore use the techniques of ^omb (1976) and S cargle 



(1982) to calculate the PSD, and we only consider frequen- 
cies / Ss 10~ 4 Hz (evenly sampled in intervals of the inverse of 
the observation duration). Higher time resolution lightcurves 
showed excess power on the w 5 ks orbital time scale. The re- 
sults, binned over the greater of four contiguous frequency bins 
or logarithmically over / — » 1.15 /, are presented in Fig. |[ 
Here we use a one-sided normalization where integrating over 
positive frequencies yields the total mean square variability rel- 
ative to the squared mean for the particular lightcurve being 
analyzed. The ASCA lightcurve shows 28% rms variability, 
whereas the RXTE lightcurve shows 16% rms variability. Both 



PSD have comparable shapes, i.e., flat from ss 10~ 6 -10~ 5 Hz 
and slightly steeper than /"' at higher frequencies. 

We have fit a broken power law to the data, assuming er- 
ror bars equal to the average PSD value divided by the square 
root of the number of frequency bins averaged over. Although 
such error estimates are only valid for averages made from 
independent frequency bins (which is not strictly true for the 
Lomb-Scargle periodogram; Scargle 1982), this should provide 
a rough estimate, especially as the lightcurves are nearly evenly 
sampled. The best fit PSD slopes are -1.3 ±0.2 for the PCA 
and —1.6 ±0.3 for ASCA (errors are A\ 2 = 2.7). The break 
frequencies are found to be (8 ± 3) x 10~ 6 Hz for the PCA, and 
(1.5 ±0.5) x 10" 5 Hz for ASCA. These break frequencies are 
consistent with the value found by McHardy, Papadakis & Utt- 
ley (1998). 

At frequencies > 5 x 10~ 4 Hz there are enough contiguous 
data segments to be able to calculate t he PSDs using stan- 
dard FFT methods ( Nowak et al. 1999 , and extensive refer- 
ences therein). In Fig. g we show the calculated, background 
and noise-subtracted, PSD for the 1.8-3.6keV and 8-15 keV 
RXTE bands. Below w 2 x 10~ 3 Hz, where signal to noise 
is greatest, both PSDs are consistent with being oc f~ 2±03 
and having 6% rms variability. From the generated back- 
ground lightcurves, we estimate that background fluctuations 
contribute at most 1.5% and 2% rms variability, respectively, 
to these PSDs at < 2 x 10" 3 Hz. As the PSDs of the back- 
ground lightcurves tend to be slightly steeper than /" 2 , there 
may be some trend for background fluctuation s to steepen the 
observed high-frequency PSDs, and, in fact, Yaqoob (1997) 
find a slightly flatter (ex /"' 4 ) high-frequency PSD for MCG-6- 
30-15. 

3. TIME DELAYS 

In order to examine time delays between the low energy 
ASCA lightcurve and the high energy RXTE lightcurve, we 
use the Z-transformed Discrete Cross-Correlation Function 



(ZDCF) o f [Alexander (1997)1 whi ch is based upon the DCF 
method of Edelson & Krolik (1988). Auto-correlation functions 



can also be computed by this procedure, and we note that PSD 
derived from autocorrelations calculated via the ZDCF yield 
identical results to those presented in Fig. We use the Monte 
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Fig. 2.— PSD constructed from ASCA 0.5-2keV (filled diamonds), RXTE 1.8-3.6 keV (open diamonds), and RXTE 8-15keV (open squares) 
lightcurves. Dashed lines show the broken power law fits discussed in the text and an f 2 power law. Background and noise have been subtracted 
from the high-frequency 1.8-3.6 keV and 8-15 keV RXTEPSD, with the residual noise levels shown as a dotted and solid line, respectively. 



Carlo methods described by Peterson et al. (1998) to assess the 
significance of any uncertainties due to flux measurement errors 
as well as uneven or incomplete sampling. 

We have previously applied these methods to a simultane- 
ous Extreme Ultravio let Explorer (EUVE )/ASCA/RXTE obser- 
vation of NGC 5548 ( |Chiang et al. 2000| ), where we found ev- 
idence for the low energy ASCA band leading the high energy 
RXTE band by 5 ks. The results for the MCG-6-30-15 data are 
shown in Fig. g, for which we have used lightcurves binned 
at 512 s resolution. A positive delay indicates that the RXTE 
light curve lags the ASCA light curve. Due to the ambiguities 
associated with interpreting cross-correlation results, particu- 
larly for such small delays, we considered three different mea- 
sures of the "lag". For each Monte Carlo trial, we estimate the 
characteristic lag by (1) fitting a parabola to the ZDCF values 
to find the location of the peak, (2) computing the centroid of 
the ZDCF over positive values bracketing the maximum value, 
and (3) using the location of the actual maximum value of the 
ZDCF. In all three cases, our simulations yield evidence for a 
positive lag at various degrees of significance: ret = 0.9 ±0.7ks, 
Tcentroid = 0.9 ± 1.3 ks, T max = 0.4+^ (90% C.L.). Although the 
fitted peak estimate is consistent with a positive lag, both the 
centroid and ZDCF-maximum estimates are formally consis- 
tent with zero giving an upper limit of r < 2ks. 

At high Fourier frequency, we have used the 1.8-3.6keV 
and 8-15 keV RXTE lightcurves discussed above to search for 
frequency-dependent time lags using standard FFT techniques. 
(A complete discussion of su ch methods, includ ing calculation 
of error bars, is presented in Nowak et al. 1999 and references 
therein). No significant time delays were found, and the 1-cr 
upper limits were 50-100 s in the 5 x 10~ 4 -2 x 10~ 3 Hz range. 
Note that there were an insufficient number of uninterrupted, 
strictly simultaneous ASCAIRXTE lightcurves to allow calcula- 
tion of their relative time delays via direct FFT methods. 

4. DISCUSSION 

The MCG-6-30-15 PSD breaks to being oc /"' at « 10" 5 Hz, 
and then breaks to being approximately oc f~ 2 between 10 -4 - 



10~ 3 Hz. This is to be compared to the Cyg X-l PSD which 
has a comparable rms amplitude and shape, and has a set of 
PSD breaks at frequ encies between 0.03-0.3 Hz and 1-10 Hz 
(Nowak et al. 1999, and references therein). The black hole 



mass in Cyg X-l is estimated to be IOMq (H errero et al. 
1995); therefore, if these break frequencies scale with mass, 
then the central black hole mass of MCG-6-30-15 could be 
as low as 10 6 M Q . NGC 5548, w hich is believed to h ave a 
central black hole mass of 10 8 Mm (fc>one & Krolik 1996|; C hi- 
ang & Murray 1996; peterson & Wandel 1999 ), shows a simi- 
lar PSD with break fr equencies at sa 6 x 10~ 8 Hz and between 
3 x 10" 7 -3 x 10" 6 Hz ( |Chiang et al. 2000| ). Again, if the break 
frequencies scale with mass, then the central black hole mass 
for MCG-6-30-15 could be several orders of magnitude lower 
than that for NGC 5548. 

To date, the most car efully studied X-ray PSD f or any AGN 
is that for NGC 3516 (JEdelson & Nandra 1999| ). (See also 
McHardy, Papadakis & Uttley 1998 who present preliminary 
results for a number of AGN, including MCG-6-30-15.) For 
NGC 3516 the PSD was seen to break from nearly flat at 
^3x 10~ 7 Hz, and then gradually steepen into an /~ L74 power 
law up to frequencies as high as « 10~ 3 Hz. NGC 3516 is 
seen to be intermediate between NGC 5548 and MCG-6-30- 
15. Based upon these measur ements (and upon other fa ctors, 



argue 
i.e., interme- 
6-30-15 discussed 



such as the source luminosity), Edelson & Nandra (1999) 
for a black hole mass in the range of 1O'-1O S M0 
diate to the masses of NGC 5548 and MCG 
above. 

The RXTE/ ASCA lag upper limit for MCG-6-30-15 is also 
i ntermediate between t h e observed X-ray l ags for Cyg X-l 
([ Miyamoto et al. 1988| ; |Nowak et al. 1999| ) and NGC 5548 
( Ichiang et al. 200Cp . The time lag observed in NGC 5548, ef- 
fectively measured on the f~ 2 portion of its PSD, is 5 ks. Time 
lags on the fl at portion of the NG C 5548 PSD could be consider- 
ably longer ( JChiang et al. 2000| ). Near the PSD break from flat 
to f~ , the X-ray time lags in Cyg X-l are « 0. 1 s, whilst on the 
f~ 2 portion of the PSD the X-ray time lags are 10~ 3 -10~ 2 s. The 
MCG-6-30-15 time lags may cover a similar dynamic range 
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Fig. 3. — Z-transformed discrete correlation function (ZDCF) of the SIS, 0.5-2 keV lightcurve relative to the 8-15 keV RXTE-PCA lightcurve 
(512 s time bins). The solid line is a parabola fit to the ZDCF values. Monte Carlo simulations of these data yield an upper-limit on the RXTE time 
lag of 2.2 ks (90% Confidence Level). 



from < 2ks (overall lag) to < 100 s (high frequency lag 2 ). 

If the characteristic variability and lag times are indicative 
of mass, then a mass as low as 10 6 M Q may be required for 
the central black hole of MCG-6-3 0-15. Assuming abolo met- 



ric luminosity of 4 x 10 43 ergs s ' (|Reynolds et al. 1997), this 



would imply that MCG-6-30-15 is emitting at 30% of its Ed- 
dington rate, which is large but still plausible. A relatively low 
central black hole ma ss would make the larg e amplitude, rapid 
variability reported by Reynolds et al. ( 1 995" much easier to un- 
derstand, whereas a mass as large as the 2 x 1O 8 M , the upper 



end discussed by [wasawa et al. (1999), seems very unlikely. 



The Jj 100 s lags seen at high frequency then likely provide an 
upper li mit to the Compton diffusion time scale (see the discus- 
sion of Nowak et al. 1999). These time scales are problematic 



for future hopes of simultaneously temporally and spectrally re- 
solving the iron Ka line in this system with Constellation-X as 
it will require w 1 ks integration times to stu dy the line profile 
of MCG-6-30-15 flYoung ; & Reynolds 2000| ). 

A 30% Eddington luminosity indicates that it is worthwhile 
to search for ultra-low frequency (/ ^ 10 Hz) variability in 
excess of « 3% rms, that is, above an extrapolation of the flat 
part of the ASCA and RXTE-PCA PSDs. 'Very high state', i.e., 
JcJ O(30% LeddX galactic black hole candidates can show PSDs 
that are flat below « 10~ 2 Hz, are approximately proportional 
to f~ 2 PSD between w 10~ 2 -10 _1 Hz, are flat again between 
ps 10 _1 -1 Hz, and then break into an f~- PSD at higher fre- 
quencies. (Specifically, see Fig. 4b of Miyamoto et al. 1991, 



which shows a 'very high state' PSD of GX339-4.) The low 
frequency portion of the 'very high state' PSD has no simple 
analogy in the (usually observed) low/hard state of Cyg X-l, 



where ultra-low frequency noise is typically associated with 
dipping activity due to obscu ration by the accretion stream (A n- 
gelini, White & Stella 1994). Previous models, for example, 
have associated 'very high state' low-frequen cy variability with 
fluctuations of a viscously unstable a-disk ( [Mowak 1994 , and 
references therein). 

This highlights the major caveat that needs to be mentioned; 
we do not know the average PSD shape nor the scaling of the 
break frequencies as a function of fractional Eddington lumi- 
nosity in either galactic black hole candidates or AGN. This 
is an especially important consideration as the mass estimates 
discussed above imply a large range of Eddington luminosity 
ratios. Considering all the evidence for rapid variability and 
extremely short time lags in MCG-6-30-15 discussed above, 
however, a low mass for MCG-6-30-15 seems to us very com- 
pelling. With the advent of the X-ray Multiple Mirror (XMM) 
mission, which has large effective area and is capable of ex- 
tremely long, uninterrupted observations, these analyses will 
become more detailed for NGC 5548 and MCG-6-30-15, and 
will allow one to develop a statistical sample of numerous other 
AGN. 
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"Although the low energy band here differed from the ASCA band, if the time lags scale logarithmically with energy ( ^owak et al. 199S and references therein) 
we would have expected the lag with respect to the ASCA band to be approximately a factor of 2 greater. 
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